Large Angular Scale CMB Anisotropy Induced by Cosmic Strings 
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We simulate the anisotropy in the cosmic microwave background (CMB) induced by cosmic strings. By numerically 
evolving a network of cosmic strings we generate full-sky CMB temperature anisotropy maps. Based on 192 maps, we 
compute the anisotropy power spectrum for multipole moments I < 20. By comparing with the observed temperature 
anisotropy, we set the normalization for the cosmic string mass-per-unit-length fi, obtaining G/i/c 2 — 1.05 Iq'^o x 10 -6 , 
which is consistent with all other observational constraints on cosmic strings. We demonstrate that the anisotropy 
pattern is consistent with a Gaussian random field on large angular scales. 
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Cosmic strings are topological defects which may have 
formed in the very early universe and may be respon- 
sible for the formation of large scale structure observed 
in the Universe today Q. In order to test the hypoth- 
esis that the inhomogeneities in our universe were in- 
duced by cosmic strings one must compare observations 
of our universe with the predictions of the cosmic string 
model. This Letter presents results of detailed computa- 
tions of the large angular scale cosmic microwave back- 
ground (CMB) anisotropies induced by cosmic strings 
g. These predictions are compared to the large scale 
anisotropies observed by the COsmic Background Ex- 
plorer (COBE) satellite. Because the predicted tempera- 
ture perturbations are proportional to the dimensionless 
quantity G/i/c 2 where G is Newton's constant and c the 
speed of light, one may constrain the value of fi, the mass 
per unit length of the cosmic strings. We believe that our 
estimate of fj, is the most accurate and reliable to date. 

Our methodology for computing the large angle 
anisotropy is to simulate the evolution of random realiza- 
tions of a cosmic string network || . From these network 
simulations we construct the temperature anisotropy pat- 
tern seen by various observers within the simulation vol- 
ume. We have evolved the strings from a redshift z = 100 
to the present, in a cubical box whose side length is twice 
the Hubble radius at the end of the simulation. This large 
box assures that the anisotropy pattern is unaffected by 
the finite simulation volume. 

In order to obtain the large dynamic range required for 
these simulations we have used a new technique whereby 
the number of segments used to represent the string net- 
work decreases as the simulation proceeds. We have con- 
ducted tests of this method by comparing smaller simu- 
lations, with and without decreasing the number of seg- 
ments: the average long string energy density is unaf- 
fected; the distribution of coherent velocities (the string 
velocity averaged over a particular length scale) is pre- 
served down to scales smaller than 1/100 of the horizon 



radius; the effective mass-per-unit-length of string (the 
energy in string averaged over a particular length scale) 
is preserved down to scales smaller than 1/100 of the 
horizon radius. The decrease in the number of segments 
was regulated so that on the angular scales of interest the 
simulation provided a good representation of the cosmic 
string network. Here we are interested in comparing with 
data from the COBE Differential Microwave Radiometer 
(DMR) which measures the anisotropy convolved with 
an approximately 7° FWHM beam B . Our contact with 
this data is through COBE's predicted correlation func- 
tion at the 10° angular scale. 

The temperature patterns are computed using a dis- 
cretized version of the integral equation 



AT 

-jr{n,x ohs ) 



d A xG^(n,x obs ,x)e^(x) (1) 



where Q^u is the stress-energy tensor of the string at 
4-position x, x Q bs is the 4-position of the observer, and 
n indicates the direction of the arrival of photons on the 
celestial sphere. Our results apply only if the present spa- 
tial curvature and cosmological constant are small since 
we use a Green's function G MI/ appropriate to a matter- 
dominated Einstein-deSitter cosmology p). While the 
universe is not entirely matter-dominated at redshifts 
close to recombination (z lec ~ 1100), the approximation 
is justified because we show below that the large-scale 
anisotropies are primarily produced at z < 100. Hydro- 
dynamical effects at recombination, which only slightly 
affect the large angle anisotropy, are not included here. 
Instead we assume that that the photons are comoving 
with the dark matter at recombination. We have used the 
ansatz of local compensation j6) as the initial condition 
for the perturbations of the matter distribution. 

We have generated three independent realizations of 
the cosmic string network. For each realization we have 
computed the fractional CMB temperature perturbation, 
■^(ri,a;obs), in 6144 pixel directions, n, on the celestial 
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sphere of 64 observers distributed uniformly through- 
out each simulation box. This computational scheme 
gives ^j!r- smoothed on about the size of the pixels which 
are 3.5° or smaller 0. One of these temperature 
maps is shown in figure [j]. This temperature map has 
been smoothed with a 10° FWHM beam to permit di- 
rect comparison with the published COBE sky maps 



(http : //www. gsfc.nasa.gov/astro/cobc/dmr jmagc.html) 
Note that these published maps have structure on angu- 
lar scales smaller than 10°; this is receiver noise. 

Each temperature map may be expressed in terms of 
the scalar spherical harmonics Yt m on the sphere: 
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In this Letter we only consider coefficients ai rn with 
£ < 20. For £ > 20 the error due to finite pixel size 
and gridding effects is larger than 10%. In reference |t]] 
we consider how to correct the large £ harmonics for dis- 
cretization effects. 

For each map we construct the multipole moments 
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The monopole and dipole moments (£ = 0, 1) are dis- 
carded because they contain no useful information. The 
mean and variance of £(£ + l)Ci for 2 < I < 20 is plotted 
in figure || in units of (G^i/c 2 ) 2 . We see that the mean 
spectrum is roughly flat, with £(£ + l)d e ~ 350(G^/c 2 ) 2 
for £ < 20. 

The standard estimator for the ensemble-averaged cor- 
relation function G( 7 ) of at angle 7 may be expressed 
in terms of the Ci. We smooth the temperature pattern 
first with the average DMR beam model window function 
Gi (tabulated values are given in which is approxi- 
mately a 7° beam, and second with a 7° FWHM Gaus- 
sian window function Wg(7°), for an effective smoothing 
on angular scales amooth = 10°. 
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Here the Pe are Legendre polynomials. By neglecting 
£ > 20 we would underestimate this sum for # amooth < 
10°. We compute the correlation G(0°, 10°) at the angle 
7 = 0° from our cosmic string simulations and compare it 
with the same quantity as estimated by the COBE team 
from the COBE data. 

We expect that most of the anisotropy on a par- 
ticular angular scale is seeded by the strings near the 
time when the projection of the coherence length of the 



string network subtends that angle on the celestial sphere 
P,f0 11 l^j . Since the coherence length of the string net- 
work grows with time, we expect anisotropics on small 
angular scales to be seeded at early times, and the large 
angle anisotropics to be seeded at late times. We have 
constructed temperature maps which include only a part 
of the temporal evolution of the string network, from a 
redshift z to the present. By considering the convergence 
of G(0°, 10°) as z increases we may determine what red- 
shift range is required to accurately determine G(0°, 10°). 
In figure H we plot the average value of G(0°, 10°) in units 
of (G/i/c 2 ) 2 for 64 observers in each of three string simu- 
lations as a function of z. We see that G(0°, 10°) receives 
its dominant contribution within a redshift z ~ 20, al- 
though there is continued growth through z ~ 100. Two 
independent models of the temperature correlation func- 
tion indicate that by neglecting the contribution of cos- 
mic strings in the redshift range 100 < z < z rcc we under- 
estimate G(0°, 10°). The underestimate in the smoothed 
autocorrelation function is 18% according to the semi- 
analytic model of Bennett et al pj], and 25% according 
to the analytic model of Perivolaropoulos |0]. We use 
the latter, more conservative model to extrapolate our 
results, valid to z = 100, out to the redshift of recombi- 
nation. 

We may normalize the cosmic string mass per unit 
length, /i, by matching estimates of G(0°,10°) from 
COBE-DMR with our predictions. The COBE-DMR 
four year sky maps yield G(0°,10°) = (29 ± ImK/T] 2 
jl3| with mean temperature T = 2.728 ± 0.002 K @. 
Our simulations indicate that at z = 100, G(0°,10°) = 
82 ± 19(G/i/c 2 ) 2 where the ±19 gives the cosmic vari- 
ance between the different observers in all three of our 
simulations. The analytic model |llj] predicts that at 
z roc , G(0°,10°) = 103±24±20(G/x/c 2 ) 2 , where the ±24 
is the new cosmic variance. The ±20 is a conservative 
systematic error composed of a ~ 10% uncertainty due 
to the simulation technique of reducing the number of 
string segments, ~ 7% due to the difference in the two 
models used for the extrapolation out to z rec , and ~ 5% 
due to the discretization of the celestial sphere. (These 
errors will be discussed in more detail in 0.) Note that 
our result is not strongly dependent on the extrapolation 
out to z rec , which makes only a small correction, lower- 
ing the normalization of fi by 10%, in comparison to the 
quoted uncertainties. Hence, adding these errors linearly, 
normalization to COBE yields 
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for the cosmic string mass per unit length. 

Because the spatial distribution of the cosmic string 
network is not described by Gaussian random variables, 
we expect that the anisotropy pattern generated in the 
cosmic string scenario will be non-Gaussian at some level. 
At very small angular scales, the sharp temperature dis- 
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continuities across strings guarantee non-Gaussian fea- 
tures. One might hope to find a non-Gaussian signa- 
ture to distinguish cosmic string models from inflationary 
models, for which the anisotropy patterns are expected 
to be very Gaussian. On the large DMR angular scales 
we are studying, however, we will see that many differ- 
ent strings contribute significantly to each resolution el- 
ement of the temperature pattern, so the conditions for 
the central limit theorem are very well satisfied and the 
temperature pattern is very close to Gaussian. 

We have looked for non-Gaussianity in the distribution 
function of the temperature anisotropy after smoothing 
our maps with the average DMR beam (s) model win- 
dow function, an approximately 7° Gaussian beam. The 
distribution function after combining all of our maps is 
shown in figure [|. We see that the distribution, on an- 
gular scales accessible to DMR, is very close to a Normal 
distribution. The pixel temperature distribution for any 
such single observer will not appear as smooth, just as for 
a limited sample drawn from a true Normal distribution. 

It has been suggested that the distribution of tem- 
perature differences is a better test of non-Gaussianity 
than the temperature distribution |I^ , |l6| ] . In figure || we 
plot the distribution of the differences in temperature of 
nearby pixels. We see that for temperature differences 
on angular separations greater than the COBE DMR 
7° angular resolution scale, the distribution of temper- 
ature differences is again very close to a Normal distri- 
bution. With finer angular resolution than that probed 
by DMR, as may be possible with the MAPS or CO- 
BRAS/SAMBAS detectors, the inherently non-Gaussian 
character of the temperature anisotropy due to cosmic 
strings may be observed. We shall further examine this 
and other aspects of the map statistics in reference |?J . 

In this Letter we have presented the first computa- 
tion of large angle CMB anisotropy from cosmic strings 
which has included all of the relevant physics. Our 
normalization of /i (|J) is consistent with most previ- 
ous work. Existing studies of the expected large-scale 
CMB anisotropics find Gfi/c 2 = 1.5(±0.5) x 1(T 6 ffij 
and 1.7(±0.7) x 1CT 6 |J (also see jjj). Note that our 
results do not appear consistent with Coulson et al |jl6| , 
who obtain the higher normalization G\ijc? = 2 x 10~ 6 
(they quote no uncertainties). They use the lattice-based 
Smith- Vilenkin evolution algorithm in Minkowski space- 
time to produce 18 realizations of a 30° square patch 
temperature field. Our results should be more accurate 
because: we simulate the string motion in an expand- 
ing universe; our simulation allows smooth variation of 
quantities such as the string's velocity; our 192 full-sky 
realizations are comparable to over 8800 30° patches and 
thus have better statistics. Our normalization of is also 
compatible with other observational constraints on cos- 
mic strings. The bound on gravitational radiation due 
to pulsar timing residuals and primordial nucleosynthe- 
sis gives Gfi/c 2 < 5.4(±1.1) x 10~ 6 Q. The bound due 



to cosmic rays emitted by evaporating primordial black 
holes formed from collapsed cosmic string loops gives 
G/i/c 2 < 3.1(±0.7) x 10- 6 CMB fluctuations on an- 
gular scales below 10 arc-seconds gives G^/c 2 < 2 x 10 -6 
at the 95% confidence level, assuming no reionization 
p9[ . The normalization presented in this Letter should 
allow a more direct confrontation of the cosmic string 
model with observations of large scale density inhomo- 
geneities. 
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FIGURES 

FIG. 1. The temperature anisotropy map for one of 192 re- 
alizations of the cosmic string induced CMB anisotropy pat- 
terns which we have generated for this Letter. This is an equal 
area, all-sky map using a Hammer- Aitoff projection. The map 
has been smoothed with a 10° FWHM Gaussian beam. 

FIG. 2. Plotted is the £(£+l)Ci angular power spectrum for 
the anisotropies from our realizations. The central dot gives 
the mean for all 192 observers while the symmetric error bars 
give the rms variation between different observers. 

FIG. 3. Plotted is the rms anisotropy after smoothing with 
an effective 10° smoothing beam versus the maximal redshift, 
2, for which strings have been included. The different curves 
give the mean value of rms anisotropy averaged over the 64 
observers in each of the three cosmic string simulations. Local 
compensation has been applied at each starting time. 

FIG. 4. Plotted is the frequency of temperature anisotropy 
for all pixels in each of our 192 maps after smoothing with a 
model of the average DMR beam, an approximately 7° Gaus- 
sian beam. The solid line gives a Normal distribution with 
zero mean and variance determined by the pixel distribution 
amplitude. 

FIG. 5. Plotted is the frequency of temperature differences 
for all pairs of pixels at a given angular separation in each 
of our 192 maps with no smoothing. The lines show a Nor- 
mal distribution with variance determined by the pixel differ- 
ence distribution amplitudes. The three pairs of curves with 
increasing variance are for separations spanning the ranges 
10° ± 1°, 20° ± 1°, and 40° ± 1°. 
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